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In this paper we study the effect of substituting R (rare-earth ion) by non-magnetic ions in the 
spin-1 chain material R2BaNi05. Using a strong-couphng expansion and numerical density matrix 
renormalization group calculations, we show that spin- wave bound states are formed at the impurity 
site. Experimental consequences of the bound states are pointed out. 
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The antiferromagnetic Heisenberg spin chains has been 
a subject of immense interest since Haldane pointed out 
that the low energy physics of integer and half inte- 
ger spin chains are fundamentally different [Q. With 
advance in both experiments and numerical tech- 
niques [Q, the prediction by Haldane is now generally 
accepted. More recently, a series of experiments on the 
family of quasi-one-dimensional materials with a gen- 
eral formula R2BaNi05 p[]5|-[ri[| where R is one of mag- 
netic rare-earth elements substituting fully or partially 
Y have made it possible to study the effect of a stag- 
gered magnetic field on a S* = 1 spin chain. The com- 
pound YzBaNiO^ is believed to be a good realization of 
5=1 Haldane gap system The S — 1 spins are 
coming from the Ni^^ ions. Other members have ad- 
ditional magnetic ions R^^{R ^ Y) positioned between 
two neighboring Ni chains. The coupling between these 
ions and Ni"^^ ions are weak. Nevertheless, these mag- 
netic ions are AF ordered below certain Neel temperature 
T/v, and the ordering affects the Ni chains by providing 
them an effective background staggered magnetic field. It 
is expected that the three modes of the Haldane triplet 
will be splitted by the staggered field at T < T/v, with 
the longitudinal mode energy increases roughly twice as 
fast as the transverse ones with decreasing temperature 
below Tat [p^[-p^ . These effects were indeed observed in 
neutron scattering experiments [p|-pd]|. 

In this paper we shall study the effect of replacing R 
ions by non-magnetic impurities. In the case of pure 
5 = 1 Haldane spin chains, the effect of replacing a spin 
by non-magnetic ion has received much attention because 
of the theoretical prediction of 5 = 1/2 excitations local- 
ized at the ends of broken 5=1 spin chains ||l|l. The 
existence of end excitations with fractional spin magni- 
tudes are believed to be a general property of Heisenberg 
spin chains and is a consequence of the topological prop- 
erties of these systems |l7]. However, the effect of sub- 
stituting R by non-magnetic impurities is very different. 
As far as the Ni^^ chain is concerned, the replacement of 
an ion by an non-magnetic impurity just alters the 
local effective (staggered) magnetic field seen by the spin 



chain, and the topological character of the spin chain is 
unaffected. Therefore, topological 5 = 1/2 end excita- 
tions are not expected to be induced by the i?-ion sub- 
stitutions. Nevertheless, we shall show in the following 
that although topological end excitations are not induced 
by _R-ion substitutions, 5 = 1 magnon bound states are 
formed around the impurity and give rise to observable 
effects on the system. 

Our starting point is the effective Hamiltonian for the 
Ni chain of the R2BaNiOz family at T < T^, 



Hbuik = J ^ Si.Si+i + /i^(— l)*5f 



(1) 



where J is the exchange constant, h — gSfisHs where 
Hs ^ Ms is the amplitude of the staggered field and 
AIs 7^ at T < T/v is the sublattice magnetization from 
the R ions. Replacing an R ion by non-magnetic impurity 
at site 1 = generates an effective local magnetic field h' . 
We shall assume a corresponding impurity Hamiltonian 

H'irn = h 5g , 

in the following, where h' — ~h for non-magnetic impu- 
rity. We shall consider h' as an arbitrary parameter, and 
study it's effect on the spin chain in the following. The 
Hamiltonian H — H^uik + Him will be studied using two 
different approaches. In the first approach, we replace 
H by the corresponding lattice non-linear-cr-model (or 
coupled-rotor model) |T^], and shall study the problem 
in the strong-coupling limit when the Haldane gap mn is 
much larger than the coupling between rotors. This ap- 
proach provides a rather simple physical picture on the 
effects of staggered magnetic field and impurity on the 
Haldane chain. The results obtained from the strong- 
coupling theory is compared with results obtained di- 
rectly from diagonalizing H with the numerical density 
matrix renormalization group (DMRG) method. Com- 
bining the results from the two methods we believe that 
our study provides a clear picture on the effect of substi- 
tuting R by non-magnetic impurities in the R2BaNiO^ 
compound. First we consider the coupled-rotor model. 
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In this treatment the low energy physics of an integer 
Heisenberg spin chain in staggered magnetic field h is 
described by the Hamiltonian 



(2) 



where Hi is a unit vector located at the link between 
site I — 1 and site i, and Li is the angular momentum 
operator conjugate to n^. A and J' are parameters to 
be determined. In the strong-coupling expansion, we as- 
sume A » J' and treat J' as a small perturbation. 
In particular, to study the one-magnon excitations, we 
shall only keep the Z = 0, 1 states (Z's are eigenstates 
of angular momentum operator L.L) and shall diagonal- 
ize the Hamiltonian -ff^ot approximately in this subspace. 
We shall neglect renormalization of ground state by two 
magnon processes in our theory. First we consider the 
limit when h = 0. 

After some straightforward algebra, we obtain the ef- 
fective Hamiltonian in the Z = 0, 1 subspace, 



H. 



(0) 
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(3) 



where we have set h = \. l{l -I- 1) = 0, 2 are the eigenval- 
ues of L.L in the Z = 0, 1 subspace. m is the eigenvalue of 
Iz, a^{am) are raising and and lowering operators with 
matrix elements < lm|a^|00 >=< 00|a„i|lm >= 1, and 
with all other matrix elements equal to zero, llm > are 
the angular momentum eigenstates. The one-magnon 
eigenstates of (||) are triplets (Z = 1, m = 0, ±1) with 
energy spectrum Ck = 2(A — ( J'/3)cos(fc)), and with Hal- 
dane gap ttih — 2(A — J' /3). 

We now consider the staggered magnetic field h. First 
we consider the effects of an effective staggered magnetic 
field heff on a single rotor. The term heffn^ couples the 
|00 > and |10 > states and the new eigenstates of a single 
rotor in the 1 — 0,1 subspace are 

\G> = /3|1,0 > +a|0,0 >, |0 >= -a|l,0 > -H/3|0,0 >, 
|1 > = |1,1 >, I -1 >= >, (4) 

where |G > is the new ground state, and |m > (m = 
0, ±1) are the new excited states with Iz = m. The coef- 
ficients a and (3 are given by 

«'-^(l + 7(/^e//)-^), fi'^lil-liheffr'), (5) 



— ) and with correspond- 



h 

where -/{heff) = ^{l + f^., 
ing eigen-energies Eq ^ A(l - 7(^e//)), Eq - Eq = 
2Aj{heff), and E±i - Eg ^ A(l + j{heff))- Notice 
that the energies of the longitudinal (m = 0) and trans- 
verse (to = ±1) magnon modes are splitted by heff. 
The staggered magnetic field also induced a non-zero 



ground state expectation value of staggered magnetiza- 
tion, < Hz >=< G\nz\G >= -^OiP- To proceed further 
we replace the rotor Hamiltonian (||) by a mean-field 
rotor Hamiltonian, 



J'^n-.n-+i - heff 



(6) 



where n; Uz > +ni and heff — h + 2 J' < Hz >. 
Treating the J'n'^.n'i^i term as perturbation in the I = 
0, 1 subspace as before, we obtain after some straightfor- 
ward algebra the mean-field equation 



1 4J' 1 



(7) 



which determines heff for given h, and with correspond- 
ing eigenvalue spectrum for the one magnon states, 

2 /' 

e±i(fc) - A(l + ^(heff)) - —a'cosik), (8) 
2 /' 

eo{k) = 2A7(/ie//) -—(«'' + P^)cosik), 

for the transverse (m = ±1) and longitudinal (to = 0) 
magnons, respectively. At small h, it is easy to show from 
(0) that heff — h/{l — 2J'/3A), and the system is sta- 
ble only when 3 A > 2 J', implying that our approximate 
mean-field treatment is only valid when mn > 2J'/3. It 
is also easy to show from (^ and (||) that the Haldane 
gaps m^^\m = 0, ±1) increase as a function of h, with 



increase in Haldane gap Am^-* ~ 2Ato^^'', in qualita- 
tive agreement with experiment and numerical DMRG 
results. 

The impurity Hamiltonian Him generate an extra term 
in the rotor Hamiltonian, 



,(±1) 



H. 



(im) _ jh' - h) 



{n'^ + L'„+nl~Ll). 



(9) 



As a result, the effective staggered magnetic field act- 
ing on the spin chain becomes position dependent, with 
h{i) = /z for i ^ 0, 1, and h{i) = {h + h')/2 for i = 0,1. 
Applying mean-field analysis as before, we obtain new 
mean-field equations for every site i in the rotor model, 

J'(< n^+i > + < >) - heffii) ~ h{i), (10) 

,heff{i). 



> ^l{heff{t)){- 
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which can be solved numerically. In addition, the 
( ^ 2^ )L^ terms in ([l^) gives rise to shift in the ener- 
gies of the |m = ±1 > states in the local rotors at sites 
1 = 0,1 with E±^ -Eg = A(1 + ^{heffm T 

The effects of impurity magnetic field on the one 
magnon spectrum can be understood quite easily. The 
mean-field equations (|l^) produce effective staggered 
magnetic field heff{iYs smaller than h (assuming h' ~ 



2 



—h) within a region of distance ^ around the impurity 
perturbed sites z = 0, 1, where ^ is the spin-spin correla- 
tion length. As a result, the local magnitude of Haldane 
gap is reduced (see Eq. (||)) in both the longitudinal and 
transverse channels and induces an attractive potential 
well of depth ~ Am^''' and range ^ ^ for the one magnon 
states. As a result bound states of magnon are formed 
around the impurity. The effective potential well formed 
around the impurity is stronger in the longitudinal chan- 
nel than the transverse channel as the effect of /le// on 
the Haldane gap is about two times larger in the longi- 
tudinal channel. Moreover, the energies of the m = — 1 
magnon states become lower than the m = 1 states be- 
cause of the additional — ^^ )-^^ coupling in -ff^oT^- 
Fig. 1^ we show the low energy spectrum obtained in our 
theory with A — 1, J' = 1.47, h = 0.2 as a function of h'. 
Notice that in the presence of single impurity, the trans- 
lational symmetry of the Hamiltonian is broken and the 
only symmetry left is the reflection symmetry upon the 
impurity site. As a result, the low-energy one-magnon 
spectrum of the system can be divided into six sectors la- 
beled by (^f^^r"*^ = 0+, 0", 1+, 1", (-1)+, (-1)", re- 
spectively where S^g^{= m) is the z component of total 
spin, and -I- or — denotes the parity of the wavefunctions. 
We see from the energy spectrum that spinwave bound 
state develops when h' < 0.2, and there is roughly one 
bound state (except the —1+ magnons at intermediate 
value of h') per sector. 

To test our theory and to get more quantitative pre- 
dictions of the impurity effect we also carry out ex- 
plicit DMRG calculations for the low energy spectrum 
of Hamiltonian H = Htuik + Him ■ We use the standard 
DMRG algorithm |l^j2^ by targeting four states and 
keeping 400 optimal states, and have studied chains with 
length up to L = 100. We have fixed h — 0.2( J) in our 
calculation, which seems to be an appropriate value for 
the R2BaNiOc, chains jl^]. The value of h' is varied to 
study the effect of impurity on the low energy spectrum. 
The numerical results of the low energy magnon spec- 
trum in the six (-S'^qj)'"""***' sectors are shown in Fig. |^. 
Our results show that in the four sectors l"*", 1~, (— 1)^, 
(—1)+ corresponding to the transverse modes, one bound 
state is induced in each sector when h' < 0.2, in agree- 
ment with our mean-field theory. The behavior of the 
bound state energy as a function of h' is also qualita- 
tively similar to those obtained in mean-field theory. In 
the 0"*" and 0~ sectors which correspond to longitudi- 
nal modes, the situation is more complicated. We find 
that more than one bound states emerge inside the Hal- 
dane gap in contrast to mean-field theory. Notice that 
although the number of bound states differs, the quali- 
tative behavior of the lowest bound state energies as a 
function of h' in the 0"*" and 0~ sectors are qualitatively 
similar to mean-field results. 

The reason why there exists only one bound state in 
mean-field theory can be understood quite easily. Our 
mean-field theory is valid only when mn > 2 J'/3, corre- 



sponding to correlation length ^ ~ 1 lattice site whereas 
^ ~ 6 lattice sites in real spin one chains. The size of 
the effective attractive potential well is of order ~ 2 lat- 
tice sites in our theory, which is much smaller than the 
size of the effective attractive potential well (^ 12 lattice 
sites) in real spin one chains. It is therefore not surprising 
that the number of bound states in real spin one chains 
is larger than one per sector but there exists only one 
bound state per sector in our mean-field theory. 

To determine the number of the bound states in the 
longitudinal sectors, we carried out a more detailed cal- 
culation for 0"*" sector by targeting ten states and keeping 
800 optimal states for h' = —1.0 and h = 0.2, the chain 
length dependence of these states are shown in Fig. ||. 
We see from the figure that as the length of the chain 
increases, more in-gap bound states are formed. In fact, 
it is difficult to determine from our numerical result the 
exact number of bound states in this sector. 

In the physical regime h' = —0.2 which corresponds to 
non-magnetic impurity, we find that bound states exist 
in the 0+,0~,l"'" and (—1)''' sectors. The bound state 
energies in the 0"*" and (-1)'*' sectors are larger and are 
roughly 0.1 x mn, where rriH are the corresponding Hal- 
dane gap energies, whereas the bound state energies in 
the 0~ and I''" sectors are much smaller (~ 0.02 x mn- 
The existence of bound states at h' = —0.2 can be seen 
more clearly from the magnetization density difference 
< Sf >i.e. — < Si >g.s. as a function of position i, 
which measures the magnetization carried by the lowest 
excited state at each sector. The results are shown in 
Fig. 1^. We can see clearly from the figures that the first 
excited states of 0+, 0~, l"*" and (—1)''' sectors are bound 
magnon states whereas the first excited states of the 1~ 
and (— 1)~ sectors are bulk states. 

The bound states can be observed in neutron scattering 
experiments on high quality crystals as an effective reduc- 
tion of Haldane gap by ~ 10 percent at T < T/v when 
non-magnetic impurities are introduced. They also show 
up as extra spectral weights with energy O.dmn in 
dynamic structure factor S{q^ w) at wave vectors q away 
from TT. 

In summary, we have studied the effect of substitut- 
ing R (rare-earth ion) by non-magnetic ions in the spin-1 
chain material R2BaNiO^. We find in numerical DMRG 
calculation that magnon bound states appear at the im- 
purity site. The bound states can be qualitatively under- 
stood by a mean-field strong-coupling expansion treat- 
ment of the spin chain and can be observed in neutron 
scattering experiments. 
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FIG. 1. The low energy spectrum in the presence of single 
impurity as a function of h' obtained in our mean-field theory 
with A = 1, J' = 1.47 and h = 0.2. 
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FIG. 2. The the low energy excitation states of the six 
sectors for staggered magnetic field h=0.2 and chain length 
L=100, the dotted line is the corresponding value of the Hal- 
dane gap for h=0.2 without impurity. 
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FIG. 3. The chain length dependence of the lowest nine 
energy levels for 0"*" sector for h=0.2 and h'=-1.0, the dotted 
line is the gap value at the thermodynamical limit. 
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FIG. 4. The magnetization density difference of the 

lowest excitation state of each sector to the ground state, 
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< SI >g.s. at the physicsl regime h' 



-0.2. 
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